#!/usr/bin/env python3
# -*- coding: utf-8 -*-

import matplotlib.pyplot as plt
plt.clf()

VE=10.0E-3 #volume à l equivalence en L
Cacide=0.10 #concentration en quantite de matiere du reactif titrant en mol/L
nOH_i=Cacide*VE #quantite de matiere initiale du reactif titre en mol
V=[i*VE/10.0 for i in range(25)] #crée 25 valeurs de volume de solution titrante ajoute
#valeurs rangees dans une liste
nOH=[]
nH3O=[]
nNa=[]
nCl=[]

for v in V :
    if v<VE :
        nOH.append(nOH_i-Cacide*v)
        nH3O.append(0)
        nNa.append(nOH_i)
        nCl.append(Cacide*v)
    
    else :
        nOH.append(0)
        nH3O.append(Cacide*(v-VE))
        nNa.append(nOH_i)
        nCl.append(Cacide*v)

plt.plot(V, nOH, 'b+-', label="OH-")
plt.plot(V, nH3O, 'go-', label="H3O+")
plt.plot(V, nNa, 'rx', label="Na+")
plt.plot(V, nCl, 'cs', label="Cl-")
plt.xlabel("Volume de solution titrante en L")
plt.ylabel("Quantite de matiere en mol")
plt.legend()
plt.show()